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ABSTRACT 

Context. The penumbra of sunspots shows a complex magnetic field topology whose three-dimensional organization is still under 
debate after more than a century of investigation. 

Aims. I derived a geometrical model of the penumbral magnetic field topology from an uncombed inversion setup that aimed at 
reproducing the Net Circular Polarization (NCP) of simultaneous spectra in near-infrared (IR; 1.56 /an) and visible (VIS; 630 nm) 
spectral lines. 

Methods. I inverted the co- spatial spectra of five photospheric lines with a model that mimicked vertically interlaced magnetic fields 
with two distinct components, labeled background field and flow channels because of their characteristic properties (flow velocity, 
field inclination). The flow channels were modeled as a perturbation of the constant background field with a Gaussian shape using the 
SIRGAUS code (Bellot Rubio 2003). The location and extension of the Gaussian perturbation in the optical depth scale retrieved by 
the inversion code was then converted to a geometrical height scale. With the estimate on the geometrical size of the flow channels, I 
investigated the relative amount of magnetic flux in the flow channels and the background field atmosphere. 

Results. The uncombed model is able to reproduce the NCP well on the limb side of the spot and worse on the center side; the VIS 
lines are better reproduced than the near-IR lines. I find that the Evershed flow happens along nearly horizontal field lines close to the 
solar surface given by optical depth unity. The magnetic flux that is related to the flow channels makes up about 20-50% of the total 
magnetic flux in the penumbra. 

Conclusions. The gradients obtainable by a Gaussian perturbation are too small for a perfect reproduction of the NCP in the IR lines 
with their small formation height range, where a step function seems to be required. Two peculiarities of the observed NCP, a sign 
change of the NCP of the VIS lines on the center side and a ring structure around the umbra with opposite signs of the NCP in the Ti i 
line at 630.37 nm and the Fei line at 1565.2 nm deserve closer attention in future modeling attempts. The large fraction of magnetic 
flux related to the flow channel component could allow to replenish the penumbral radiative losses in the flux tube picture. 
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1. Introduction 

It has not been possible up to now to unambiguously deter- 
mine the topology of the magnetic fields inside the penum- 
bra of sunspots directly from spectroscopic or spectropolari- 
metric observations due to the many different possibilities of 
interpreting the data. The observations have provided some 
boundary condi tions, like for instance the presence of the 
Evershed flow (lEvershedl Il909l) . an almost radial orientation 
of the magnetic field lines, a radial decrease of magnetic field 
strength, and a r adial increase of the field inclination to the lo- 
cal vertical (e.g.,lLites et al. 1993; Westendorp Plaza et al. 2001; 
Borrero & Bellot R ubio 2002; Solanki 2003: Bell ot Rubio et al.1 
2004; Langhans_eiaL||2005t |BecJd|2008f Schlichenmaier 2009). 
With increasing spatial resolution, the Evershed flow could be 
definitely be shown to be related to the penumbral filaments 
(IShine et al.lll994t iTritschler et al.l 12001 lLanghans etaD 120051: 
iRimmele & Marino1l2006l) . The penumbral filaments have a very 
peculiar appearance in data of highest spati al resolution, a 
dark core flanked by two lateral brightenings (iScharmer et al.l 
2002; Sutterlin et al. 2004), with indications that the flow ve- 
locity is highest in the darkest part (IBellot Rubio et al.l [2005; 
lLanghans et al .120071: iBellot Rubio et al.l2007l) . The small scales 



of this internal structure of penumbral filaments presumably 
contributed to the fact that conflicting results on the rela- 
tion between flow veloc ity and intensity were found in sev 



eral previous studies (e.g. JDegenh ardt & Wiehr 1991 1:lTitle et al 



1993:lLites etal.1119 93: Shin e et al.lll994l: ISchlichenmaier et al 
2005; Langhans et al] l2005h . It has not been possible to 
put together all the information on the penumbra into an 
unique commonly accepted model, also thanks to another 
source of information which spectropolarimetric observations 
provide, the Net Circular Polarization (NCP), a measure of 
the asymmetry of the Stokes V polarization signal. The 
penumbra of sunspots sho ws one of the largest NCP val- 
ues o f all solar structures (lllling et al] 1 19741: lAuer & Hea slev 
19781: lHenson_^Kemrj|jl9 84: Sanchez Alme ida & Litesl 1 19921: 



Miiller et al. 



2006; Ichim oto et all l2008h . Since the NCP is 
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related to gradients of magnetic fi elds strength and velocity 
along the line of sight (LOS) (e.g., ISkumanich & Litesl [19871: 
Sanchez Almeida & Litesl! 1992), it is a crucial source of infor- 
mation on the three-dimensional organization of the penumbral 
magnetic fields. Unfortunately, gradients can come in even more 
different shapes than constant values, opening up even more pos- 
sible configurations to reproduce the observations. 

One of the suggested configurations for the three- 
dimensional organization of the penumbral magnetic fields is the 
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so-called "uncombed" penumbra (iSolanki & Montavonlll993l) . 
where less inclined fiel d lines wind around horizontal flow chan- 
nels (see for instance iBorrero et al.ll2008|). This configuration 
has been taken up for mod e ling ( Thomas & Montesinosll 1993k 



Schliche nmaier et al I [19981: ISch lichenm aier & Co llados 2002; 
Miiller et al. 2002; Borre ro et al.1 12004 [Miller et al. 2006J3 



and inversion schemes since it is one o f the few approaches 
that reproduces the observed NCP (e.g. , Martinez Pillet 2000; 
Borr ero et all2006l: lBeck 2006; Jurcak et al. 2007; Bor rero et al.1 



2007b . Another approach reproducing the observed NCP is 



the so-called micro- structured magnetic atmosphere hypothe- 
sis (MISMA; ISanchez Almeidalll998l: ISanchez Almeidal l200ll 
2005) that proposes a fine structure of the magnetic field at scales 
of a few km. On an average over the scales that are typical for 
the spatial resolution of even the most recent observations, the 
MISMA can be r epresented by two comp onents similar to the 
uncombed model (Sanchez Almeida 2005), but with additional 
properties that come from its substructure. 

Whereas the uncombed models were derived in close con- 
nection to, or in some cases, from spectropolarimetric obser- 
vations of sunspots, some other explanations for the penum- 
bral structure have been brought forward fro m a more theoret- 
ical point of view. Th omas & Weiss ( 2004) suggested the ef- 
fect of magnetic flux pumping as acting agent of the penum- 
bral fine- structure. This provides a driver for the organization of 
the penumbral magnetic fields, but gives, however, no des crip- 
tion of the organization itself. IScharmer & Spruitl (|2006) and 
Spr uit & Scharmerl (120061) proposed the existence of field-free 
gaps with convective motions below the visible surface to bal- 
ance the radiative energy losses of the penumbra, but it has not 
been shown that this model finally reproduces the o bserved NCP. 
The num erical 3-D MHD simulations described by iRempel et al.l 
(2009a b) showed the effects of magneto-convection in the sim- 
ulated sunspot, i.e., convection cells oriented and shaped by 
the direction of the magnetic field lines rather than completely 
field-free gaps, but again the cross-check of the resulting spectra 
vs. spectropolari metric observations i s still missing besides an 
initial attempt bv lBorrero et al.1 (l2010h . 

In this contribution, I investigate the NCP in simultane- 
ous spectra in the visible (VIS; 630 nm) and near-infrared (IR; 
1.56/zm) using an uncombed inversion setup that includes gra- 
dients along the LOS in field strength and velocity. The observa- 
tions are described in Sect. [2] The inversion method is explained 
in Sect. [3) its results are described in Sect.@] Section[5]discusses 
the findings, while Sect. [6] gives the conclusions. Appendix [A] 
shows the differences between initial and best-fit model atmo- 
spheres, Appendix |B] several examples of observed spectra with 
the corresponding best-fit profiles. 
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Fig. 1. Overview of the observation. Bottom row, left to right. 
continuum intensity in IR, integrated absolute Stokes V signal of 
1564.8 nm, same for 630.25 nm. Top row, left to right: line-core 
velocity of 630.37 nm, NCP of 1564.8 nm, same for 630.25 nm. 
Tick marks are in arcsec. The black dashed line denotes the lo- 
cation of a cut along the neutral line of Stokes V. The white ar- 
row points towards disc center. The black crosses and the corre- 
sponding numbers in the map of 630.37 nm denote the locations 
of the profiles shown in Figs. IB.ll to lB.3i 



calibration routines for both instruments, including t he correc- 
tion fo r the instrumental polarization of the VTT (e.g lBeck et al.l 
2005a : b). The spectra were aligned and brought to an identi- 
cal spa tial sampling as described in the Appendix of iBeck et al.l 
(|2007|). The spa tial resolution was estimated to be around 1" 
(IBeck et al.ll2007l) . In the inversion of the spectra, the transitions 
of the spectral lines and the adopt ed rest wavel engths were iden- 
tical to those given in Table 1 of lBeckl ([2008), but without the 
two weakest lines (Fe i at 1 564.74 nm and Fe i at 630.35 nm). 

The sunspot NOAA 10425 was located at an heliocentric an- 
gle of around 30°. Figure Q] shows the IR continuum intensity 
map, the integrated absolute Stokes V signal of the two more 
magnetic sensitive IR and VIS lines (1564.8 nm, 630.25 nm), the 
line-core velocity of the weak Ti i line at 630.37 nm, and the NCP 
of the lines at 630.25 nm and 1564.8 nm. The Evershed effect and 
its filamentary structure can be seen in the velocity map of the 
Tii line. This line is only slightly sensitive to magnetic fields, 
and hence allows to recover reliable line-core velocities also in 
the umbra where the other lines spli t. The two NCP map s show 
the symmetry pattern described by Miille r et all (120021 2006): 
the VIS is symmetric to the line of symmetry through sunspot 
center and disc center, the IR anti- symmetric with twice the fre- 
quency. Comparison of the NCP with the integrated V signal 
shows that the NCP is not changing its properties in the neutral 
line of Stokes V on the limb side even if the absolute V signal 
reduces strongly. 



2. Observations, data reduction and data alignment 

The observations were taken on 2003 Aug 9 with the 
POlarimetric LIttrow Spectrograph (POLIS; Bec k et al Il2005bh 
and t he Tenerife Infrared Polarimeter (TIP; lMartinez Pillet et al.l 
1 19991) at the German Vacuum Tower Telescope (VTT) in Izana, 
Tenerife, Spain. These instruments observed the Stokes vector at 
1565 nm and 630 nm, respectively. Both instruments were fed si- 
multaneously using an achromatic 50-50 beamsplitter The setup 
and the data are described in detail in Bec k et all (|2007) and 
Beck ( 2008); here I only used one of the two observations de- 
scribed in the latter. The data were reduced with the respective 

1 Thomas & Montesinos actually predates the paper of Solanki & 
Montavon. 



3. Data analysis 

The co-spatial spectra of all observed VIS and IR lines were 
first invert ed simultaneously with the standard version of the 
SIR code dRuiz Cobo & del Toro Iniestal 119921) . This analysis 
used two magnetic field components in each pixel inside the 
penumbra where the field properties were assumed to be con- 
stant with optical depth (termed "2C inversion" in the follow- 
in g). Th e inver sion setup and its results are described in detail 
m iBeckl (|2008). This inversion setup is unable to produce any 
NCP in its resulting synthetic spectra, but is, however, impor- 
tant both for comparison and as initial model for the more com- 
plex "Gaussian inversion" described below. The 2C inversion is 
not reflecting the topology of the magnetic fields in the penum- 
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Fig. 2. Example of the model atmosphere for the Gaussian in- 
version. Top row, left to right, temperature in Kelvin, electron 
pressure in dyncrn" 2 , field strength in Gauss. Bottom row. line- 
of-sight velocity in m s" 1 , field inclination to the LOS in degree, 
field azimuth in degree. Dashed: initial values of the background 
field atmosphere, solid: final best-fit background field with per- 
turbation added. The values are given vs. logarithmic optical 
depth, the final background field values can be simply read off 
from the values at log r = -4 since they are constant with depth. 



bra satisfyingly, because it corresponds to two magnetic field 
components with different field inclinations that are horizontally 
separated. To explain the observed NCP, a vertical interlacing of 
the fields has to be assumed. In the "uncombed" picture, the ge- 
ometry is commonly modeled as horizontal flow channels of a 
limited width around which the less inclined field lines wind. A 
LOS that passes through such a configuration will thus encounter 
first the less inclined field, then at some depth the flow channel, 
and eventually again the background (bg) field, depending on 
the opacity of the flow channel (fc). One way to mimic such a 
geometry in an inversion of spectra is to introduce a localized 
perturbation at some optical depth in the atmospheric stratifica- 
tion. 

To this extent, the SIRGAUS code dBellot Rubioll2003h was 
used. SIRGAUS is based on the SIR code and allows to add a 
Gaussian perturbation to otherwise constant atmosphere param- 
eters (termed "Gaussian inversion" in the following). The pertur- 
bation is specified by its location and its width in optical depth, 
which are identical in all atmospheric quantities, and a variable 
amplitude of the perturbation in each parameter. However, still 
two magnetic components are used: one component has constant 
properties with depth, and the second component consists of the 
latter plus the added perturbation. This allows to model that the 
LOS may encounter an uncombed topology only in a part of the 
pixel, which could well happen at the spatial resolution of the ob- 
servations. The respective fill factor of both components is also 
a free parameter in the inversion. Figure [21 shows how the two 
components of the Gaussian inversion look like in the best-fit so- 
lution of one set of penumbral Stokes profiles. The total number 
of free parameters is 8 for the constant background field com- 
ponent (field strength B, inclination y, azimuth 0, LOS velocity 
v, 3 nodes for the temperature stratification T, a constant micro- 
turbulent velocity v m ; c ); the Gaussian perturbation adds again 8 
more free parameters (amplitude of Gaussian in B, y, cp, v, T, v miC , 
location r center , width <x). The additional parameters shared by 
both components are an identical macro-turbulent velocity v mac , 
the fill factor of each component /, and the stray light con- 
tribution to the spectra, J3. This gives in total 19 free parame- 



ters for 329 wavelength points x 4 Stokes parameters (= 1316 
measurement values), whic h, however, are not all f ully indepen- 
dent (see the discussion in San chez Alme ida 200 5]). SIRGAUS 
has been used in a number of investigations dBellot Rubiol 
2003; Beck 2006; Jurc ak et al.ll2007l:ICabrera Solana et all2 008: 
Ishikawa et al. 20 Kj. The main difference between the standard 
version of SIR and SIRGAUS is that spectral lines cannot be 
treated as blends of each other, i.e., the mutual influence of the 
line wings on the neighboring lines is not included. For that rea- 
son, also the two weakest lines in the observed spectra (Fe i at 
1564.74 nm and Fe i at 630.35 nm) were not used for the fit since 
they are located in the wings of stronger lines. 

I found that the Gaussian inversion is very sensitive to the 
initial model atmospheres, i.e., the code fails to converge if the 
observed profiles are too far off from those resulting from the ini- 
tial model. To overcome this difficulty, I used the results of the 
2C inversion as input. The initial amplitude of the Gaussian was 
derived by subtracting the best-fit values of the atmospheric pa- 
rameters of the two components in the 2C inversion. This intro- 
duces some bias in the inversion results, since the 2C inver sion 
alread y gives a reasonable fit to the observed spectra (cf . iBeckl 
2008). However, the main properties of the field components like 
field strength and field inclination are "robust" quantities when 
using the near-IR spectral lines. These lines give rather strict lim- 
itatio ns on for instance the field strength due to their strong split- 
ting dBeck et alj|2007l Fig. 7). The Gaussian inversion also is 
aimed at retrieving the vertical organization of the field compo- 
nents rather than their basic properties. The initial location and 
width of the Gaussian perturbation were set to logr = -0.5 and 
0.5 units of log r, respectively. The same initial values were used 
for every inverted pixel in the penumbra, since these are exactly 
the quantities that the inversion should determine. 

Appendix lAl shows a comparison between the initial model 
atmospheres as derived from the 2C inversion and the final best- 
fit solutions of the Gaussian. The basic properties of the mag- 
netic vector field (field strength and orientation) changed only 
slightly. Figure IA.1I shows that the field topology corresponds 
to the predictions for an uncombed model with horizontal flow 
channels: the difference of field inclination between the two 
components, Ay, is positive throughout the spot, whereas the dif- 
ference of field azimu th AO changes s ign across the symmetry 
line of the spot (cp. iMuller et al.ll2002l their Fig. 11). 

4. Results 

4.1. Individual profiles and fit quality 

Figure [3] shows the observed and best-fit profiles of both the 2C 
and the Gaussian inversion on a cut along the neutral line of 
Stokes V (marked in Fig. Q} for a visual control of the qual- 
ity of the inversion. Each row corresponds to the spectrum of 
a single pixel. The spectra have been normalized separately to 
their respective maximal value in IQUV to improve the vis- 
ibility of spectral features. The polarization signal of the Tii 
line at 630.37 nm is close to the noise level outside the um- 
bra, it only shows up clearly in Stokes V at the locations where 
the splitting of for instance 1568.4 nm is largest (y ~ 20"). 
The complex shape of the profiles can be seen in the Stokes 
V graph (lower right). The V signals of all lines often show a 
local minimum and a maximum in the blue lobe (black/white), 
and in the red lobe either a m inimum or a maximum (see also 
iFranz & Schlichenmaierl2010h . The cut actually crosses the neu- 
tral line of the sunspot, the main polarity of the Stokes V signal 
changes at around y ~ 12". Both inversion setups reproduce 
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Fig. 3. Spectra of a spatial cut along the neutral line of Stokes V. Clockwise, starting left top: Stokes IQVU. Each subpanel shows 
the observed spectra in the top row, the best-fit profiles of the Gaussian inversion in the middle row, and those of the 2C inversion 
in the bottom row. Left to right in each subpanel: 1564.8 nm, 1565. 2nm, 630.15 nm, 630.25 nm, 630.37 nm. The x-axis (dispersion) 
is in mA, the y-axis (spatial position along the cut) in arcsec. All profiles were normalized individually to improve the visibility of 
the spectral patterns; the range is ±1. 



the observed Stokes spectra fairly well; neither differences be- 
tween observed and best-fit profiles nor between the two sets of 
best-fit profiles can be easily discerned. One case of a clear dif- 
ference can be found in the Stokes V signals of the VIS lines 
at 630.15 nm and 630.25 nm (3rd and 4th column in the lower 
right panel). In the lower part from y = 0" to 12", the observed 
V signals show a (weak) minimum (black) followed by a broad 
maximum (grey to white) in the blue lobe, and a broad minimum 
(black) in the red lobe. The 2C inversion is unable to clearly re- 
produce the broad minimum in the red lobe. Another deviation 
between observations and best-fit profiles is seen in Stokes Q of 
again the VIS lines: in the same lower part of the spectra, the 
Stokes Q of the Gaussian inversion shows strong pixel-to-pixel 
variations in, e.g., 630.15 nm with Doppler excursion to the blue 
and to the red that are absent in the observations. Appendix |B] 
contains six examples of individual observed profiles, with the 
best-fit profiles of both inversion setups overplotted. The loca- 
tion of the profiles is marked in Fig.[T] 

To quantify the differences between the two inversion setups, 
I calculated the x 2 values of the squared difference between ob- 
served (S j r0 b s ) and best- fit profiles (S j rSy nth) for both inversion 
setups by 



Xij 



hobs $ j,synth) (Ai) 



where i cycles through the spectral lines and j through the entries 
of the Stokes vector, S . 

The bottom row of Fig. shows the difference of Xy 

be- 
tween the 2C inversion and the Gaussian inversion for the spec- 
tral lines at 630.25 nm and at 1564.8 nm. Blue color indicates 
a smaller x 2 m me Gaussian inversion (=improved fit qual- 
ity), red color a larger value. For 630.25 nm, the Gaussian in- 
version improved the fit quality on the limb side of the spot 
(lower half), whereas for 1564.8 nm there are no visible changes 
there. On the center side in the innermost penumbra, Xy S ot 
worse in the Gaussian inversion (upper half). The same pro- 
files in the inner center side penumbra already show a worse 
fit in the 2C inversion as compared with profiles on the limb 
side. The reason is that the profiles on the limb side show 
clear signatures of two different magnetic field components (par- 
allel and anti-parallel to the LOS) that contribute to the ob- 
served spectra, especially in the V sp ectra from near the neu- 
tral line with more than two lobes (e.gjGrigoriev & Katzll 1 972k 
ISchlichenmaier & Colladosll2002t ISanchez Almeidall2005l) . ~The 
profile shape of the near-IR lines varies even stronger than that 
of the 630nm lines (see Fig s. IB.ll to lB.31 or Riie di et aD 



1999; 



2002; 



(1) 



del T oro Inie sta et all 120011; ISchlichen maier & Colladosl ._ 7 

Bec]d|2008). On the limb side, this signature gets lost due to the 
fact that the field lines of both magnetic components are parallel 
to the LOS. The spectra there seemingly do not provide enough 
information to fix the free parameters in the Gaussian inversion. 
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Fig. 4. Xy °f the two inversion methods. Top row: Xy of 
630.25 nm in the 2-component inversion (left), same for 
1564.8 nm (right). Middle row: same as above for the Gaussian 
inversion. Bottom row: difference of Xy between the Gaussian 
and 2C inversion. Negative values indicate a better fit quality of 
the Gaussian inversion. 



del Toro Iniesta et al.l (l201Qh in a recent publication suggest to 
use Occam's razor principle for such a case: reject the more 
complex solution (the Gaussian inversion) in exchange for the 
simpler one (2C inversion). This approach, however, is not nec- 
essarily justified even in the case that the more complex solution 
gives a worse x 2 - The 2C inversion has intrinsically a zero NCP, 
whereas the Gaussian inversion actually still reproduces the ob- 
served NCP to some degree, even at the cost of a worse x 2 - The 
more complex solution therefore should be preferred over the 
simpler one for being a more realistic model approach for the 
presumably "real" solar magnetic field topology. 

4.2. Net Circular Polarization 
The NCP was defined by 



NCP 



/ 7< 



(A)dA . 



(2) 



Note that V(A) is not normalized with the continuum inten- 
sity, I c , but the intensity at the corresponding wavelength, 1(A). 
Especially for the deep VIS lines the normalization by I c can 
change the NCP value significantly. I applied Eq. ([2]) to all spec- 
tral lines in the observations and the best-fit profiles separately, 
restricting the integration range [Aq, A\] to encompass only a sin- 
gle line each time. I remark that the NCP is no quantity, whose 
deviation between observed and synthetic spectra is minimized 
in the fit procedure. The SIRGAUS code minimizes the total 
deviation, x 2 t = YuijX 2 ^ between the observed and the syn- 
thetic profiles. If an actual misfit of the NCP value reduces x 2 v 
the code is nonetheless forced to choose the solution with the 
"worse" NCP. With Eq. ©, the derivation of the NCP is straight- 
forward, but I think that a word of caution is appropriate. The 
NCP is a favorite toy for theoretical considerations, but from an 
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Fig. 5. Comparison between the NCP in the observed (left col- 
umn) and best-fit profiles (right). Top to bottom: 630.37 nm, 
630.25 nm, 630.15 nm, 1565.2 nm, 1564.8 nm. For the fit result 
of 630.37 nm and of 1565.2 nm, the display range is half that of 
the observed NCP. 



observational point of view it is an ill-determined quantity. The 
NCP corresponds to a subtraction of, even with the integration in 
wavelength, two "small" numbers, the areas below the circular 
polarization lobes. Due to the low polarization signal level of at 
maximum about 40 % of /, the NCP value is severely affected 
by the noise level in observations. Thus, even if the NCP is one 
of the few quantities directly dependent on the vertical structure 
of the solar atmosphere, its absolute value has to be treated with 
care whenever observations are concerned. 

Figure [5] shows a direct comparison of observed and best- 
fit NCP for the whole sunspot. The display range for the best- 
fit NCP had to be halved for 1565.2nm and 630.37 nm to al- 
low seeing the spatial variation across the spot at all. The (anti- 
symmetry properties of the NCP are recovered completely only 
for the VIS lines, with positive (negative) NCP values on the 
(center) limb side. The best reproduction of the large-scale spa- 
tial pattern of the NCP is found for 630.37 nm. For 630.15 nm 
and 630.25 nm, the center side is less well reproduced. The neg- 
ative NCP values at the outer penumbral boundary in the obser- 
vations are missing in the best-fit NCP, and the values in the mid 
and inner center side penumbra are larger and strictly positive 
in the best-fit NCP. For the IR lines, the agreement between ob- 
served and best-fit NCP is not much better. The antisymmetric 
pattern of two minima and maxima on an azimuthal path can 
be seen in the best-fit NCP for 1565.2 nm and 1564.8 nm, but 
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Fig. 6. Scatterplots of NCP on the limb-side in the observed and best-fit spectra. Grey dots show all data points, black pluses the 
same with binning. The grey dashed lines are linear regression lines with the results for offset and slope as given in the lower right. 



the NCP amplitude is off by a factor of 2 for 1565.2 nm like for 
the Ti i line. Again, the agreement on the center side is gener- 
ally worse than for the limb side. An interesting feature is seen 
near the umbral-penumbra boundary. The umbra is encircled by 
a ring of negative (positive) NCP in the 1565.2 (630.37) nm line 
in the observations, which does not appear for the other spectral 
lines. The best-fit NCP actually reproduces this feature, with op- 
posite signs in the two spectral lines, even if this includes exactly 
those locations where the Xy of 1564.8 nm and 630.25 nm was 
worse. 
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To investigate how far the inversion also catched the abso- 
lute value of the NCP, I show scatterplots of observed vs. best-fit 
NCP values for all spectral lines (Fig. [6]). I restricted the area 
considered to the limb side of the spot since the 2-D maps of 
Fig. [5] already suggest that no clear correlation can be expected 
on the center side. The later detailed investigation of the inver- 
sion results also will focus on the limb side. For 630.15 nm and 
630.25 nm, the relation between observed and best-fit NCP is 
close to linear with a slope of about unity. The best-fit NCP val- 
ues are slightly overestimating the observed NCP for 630.15 nm 
(slo pe 1.21) and underestimating it for 630.25 n m (slope 0.84) 
(cp. ISanchez A lmeida 2005: lBorrero et al . 2006, Figs. 18 and 5, 
respectively). For both near-IR lines and 630.37 nm, the fit NCP 
is significantly smaller than the observed NCP, with only a weak 
trend between observed and best-fit NCP. 



Fig. 7. Radial variation of the observed NCP {black), best-fit 
NCP {green), and the LOS velocity {orange) on the limb side 
{solid lines) and the center side {dashed lines). The vertical solid 
lines denote the average inner and outer penumbral boundary. 



4.3. Radial variation of NCP and v LO s 



iTritschler et al] (120071 T07) pointed out the existence of a zero- 
crossing of the NCP of the 630.25 nm line on a radial cut on 
the center-side, where the NCP changes sign in the mid to outer 
penumbra. The maps in Fig. [5] support this finding, which also 
applies to the 630.15 nm line as well. For a direct comparison 
with their Fig. 3, I determined the radial variation of the NCP 
and the LOS velocity for the 630.25 nm line (cf. Fig. [7]). The 
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zero-crossing of the center side NCP is located at a similar ra- 
dial position of around 0.8 r/r spot as in T07, presumably to the 
similar heliocentric angles of the observations (30° and 42°, re- 
spectively). The location of the zero-crossing is co-spatial to 
the maximum of the LOS velocity on the center-side (orange 
dashed line). Figure [7] also demonstrates that the best-fit NCP of 
630.25 nm on the limb side follows closely the observed NCP 
qualitatively and, in this case, even quantitatively. 

4.4. Topology of the flow channels from the location and the 
width of the Gaussian perturbation 

The information on the location and spatial extent of the flow 
channels is contained in the central position r center and the width 
o~ of the Gaussian perturbation. In the following, I take the full 
width at half maximum (FWHM) of the perturbation as a mea- 
sure of the vertical extent of the flow channels. Figure[8]displays 
the azimuthally averaged values of location and width in optical 
depth units as function of the radial distance to the spot center, 
separately for the center and the limb side. For r/r spot < 0.6, the 
center of the perturbation always stays close to the log r = 
level. The width is so large that the lower boundary is close to 
or below the log r = level. Around r/r spot = 0.65, there is a 
distinct maximum in the location value. The maximum is more 
pronounced in the limb-side. The half-width-at-half-maximum 
(HWHM, top panel) decreases from about 2 to 1 units of optical 
depth after the maximum in the location. For r/r spot > 0.8, the 
location is roughly constant at log r = -0.5 on the limb side and 
at log r = on the center side, respectively. 

For an interpretation of the behavior one has to take into ac- 
count that the values of both parameters are given in the optical 
depth scale. If the Gaussian perturbation is located at log r = 0, it 
can be interpreted in two ways: it may indicate a low-lying flow 
channel or indicate that it is an optically thick structure. Thus, 
the isolated maximum at r/r spot =0.65 can either indicate flow 
channels sinking (location in logr changes from to +1.5) and 
immediately again rising (location in logr changes from +1.5 to 
0) in the mid penumbra, which seems improbable, or that they 
are optically thick up to a radius of r/r spot = 0.65. The follow- 
ing decrease of the location from logr = 1.5 to around -0.5 with 
increasing radial distance suggests in my opinion that the latter 
interpretation is true. As long as the LOS can not penetrate down 
to the lower boundary of the flow channel, the inversion code has 
no information on where to place it. The decrease of the location 
for larger radii could then be interpreted as the rise of a flow 
channel through the surface, where the profiles then contain the 
information necessary to place it in the atmosphere. Another ar- 
gument for this interpretation are the relative fill factors of the 
two inversion components. At r/r spot = 0.6, the fill factor of the 
stronger inclined c omponent identified with the flow ch annels 
strongly increases dBellot Rubio et al]l2004 lBeckll2008h . This 
indicates that the flow channels start to have a stronger contribu- 
tion to the observed profiles which could be due to the fact that 
they breach the surface there. It also indicates that for the inner 
penumbra the location of the Gaussian perturbation derived by 
the inversion code may be unreliable. In the inner penumbra, the 
Gaussian inversion reproduces the NCP with the discontinuity 
at the upper boundary of the flow channel, but it cannot locate it 
in optical depth because the LOS does not penetrate through the 
flow channel. 

It is worthwhile to compar e Figs. [8] with the corresponding 
Fig. 6 of iBorrero et al.l ([2006). The flow channels in the cur- 
rent investigation are globally located at lower optical depths, 
which could be the result of the initial atmosphere stratifica- 
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Fig. 8. Parameters of the Gaussian perturbation. Top: Radial vari- 
ation of the HWHM of the Gaussian perturbation in the inversion 
in units of log r. Middle: location of the center of the pertur- 
bation for the limb-side (thick line). The thin lines indicate the 
FWHM. Bottom: same as above for the center- side. 
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Fig. 9. Conversion curve between optical depth, r and geometri- 
cal height, z, from the HSRA model. +: HSRA, solid line: poly- 
nomial fit of 3th order. 



tion used for the back ground component (penumbral model 
ofldel Toro Iniesta et al.l Tl994) vs. HSRA model, respectively). 
Otherwise, the radial variation follows similar trends: the loca- 
tion is near log r = in the innermost penumbra, moves upwards 
in optical depth at around 0.7 r/r spot and shows a slight decrease 
at (beyond) the outer sunspot boundary. 
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Fig. 10. Topology of the average flow channel from the Gaussian 
inversion in geometrical height (limb side only). The inclined 
black line starting at r ~ 10 Mm gives the integrated inclination 
of the background field. The direction of the arrows gives the 
field inclination in the 2C inversion, their length is proportional 
to the velocity. Crosses mark the location of the center of the 
flow channel. 



Conversion to geometrical height I converted the optical depth 
scale to a geometrical height scale using the tabulated values of 
the Harvard Smithson ian Reference Atmosphere (HSRA) model 
(Ginge rich et al.|[l97ll) for the conversion from r into km above 
^500 = 1 (cf. Fig.©. Even if the HRSA model atmosphere has 
been derived as an average stratification of the quiet Sun devoid 
of magnetic fields, it should still be usable as first approximation. 
I assumed additionally that the iso-surface of T500 = 1 has a slope 
of 3 deg when going from the inner penumbral boundary towards 
the outer boundary corresponding to a Wilson depres sion of 380 
km in the umbra (|Schlichenmaier & Sch midt 2000). Recently, 
Puschma nn et al.l (1201 Oh have presented a more solid determina- 
tion of a geometrical height scale for inversions of spectropo- 
larimetric data that eventually could be employed to the present 
results as well. Figure [TO| shows how the location and width in 
optical depth appears in geometrical height with the assumptions 
above. For a comparison, I overplotted arrows that give the field 
inclination of the flow channel component in the 2C inversion. 
The radial variation of the location of the Gaussian perturbation 
agrees quite well with the field inclination from the 2C inversion. 

The results on the location of the perturbation seem to be un- 
reliable for r/r spot <0.65 (vertical black line in Fig. [TO]), because 
of the counterintuitive descend and rise of the location of the 
perturbation around that radius (cf. Fig. [8]). For r/r spot <0.65, 1 
thus substituted the values derived from radiall y integrat ing the 
field inclination to the surface like described in lBeckl (120081) . The 
location of the flow channel left of the black vertical line thus is 
based on the 2C inversion results. The final result of combining 
the 2C and Gaussian inversion for the whole penumbra then is 
a flow channel that ascends steeply in the inner penumbra, turns 
into a slightly elevated flow channel which is close to horizontal 
in the outer penumbra, and slightly bends down in the outermost 
penumbra. 

The conversion to the geometrical height scale also allows 
to derive the diameter of the flow channel in km. The non-linear 
relation between optical depth and geometrical height (Fig. [9]) in 
combination with a FWHM given in units of logr by the inver- 
sion makes the flow channel asymmetric around the central lo- 
cation of the perturbation in the geometrical height scale. I then 
defined the FWHM in km as the difference between the location 
of the upper and lower boundary that resulted from transforming 
T center ±HWHM (in units of logr) to the geometrical scale. The 
inversion code determined the FWHM along the LOS which can 
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Fig. 11. FWHM of the average flow channel on the limb side 
(solid). The dotted line shows the FWMH after correction for the 
field inclination to the LOS. The dashed vertical line denotes a 
FWHM of 250 km. 

cut the assumed flow channel at some angle (see Fig. [12]). The 
FWHM should thus be corrected for the field inclination to the 
LOS by a multiplication with cos (|90° - ywsjA)- Since ywsjc 
is close to 90 deg on the limb side, the effect is, however, small 
(dotted line in Fig. Qj]); I have thus not re-drawn Fig. [TO1 with 
the corrected FWHM. Within all the limitations imposed by the 
analysis and the conversion to geometrical height, the results 
indicate a rather constant diameter of the flow channel with a 
FWHM of about 250 km (cf. Fig.\TQ. 

4.5. Magnetic flux of the two components 

The relative amount of magnetic flux in the two magnetic com- 
ponents representing the flow channels and the background field 
is of interest for the question if the penumbral heating can be 
achieved by the repetitive rise or motion of hot flow channels. 
The calculation of the magnetic flux from the inversion results 
is, however, not straightforward. Figure[l2]visualizes the way the 
inversion code has constructed its synthetic spectra: the 3-D vol- 
ume element given by the spatial extent of the pixel (Ax, Ay) = 
(slit width, spatial sampling along the slit) and the optical depth 
axis is first separated in two different regions by the fill factor 
of the magnetic components, /. The fill factor is defined in the 
plane perpendicular to the LOS. Two different atmospheric strat- 
ifications along the LOS are then assumed in the two regions: in 
one case only the constant background field, in the other case 
the bg field plus the Gaussian perturbation. For the second com- 
ponent with the Gaussian perturbation, the plotted cylinder is 
actually not fully correct, since the inversion code used in princi- 
ple only a Gaussian shape along the optical depth axis, not in the 
spatial dimension. Both bg field and fc field can have an arbitrary 
orientation to the LOS and to each other; I only sketched three 
cases with a different field inclination of the fc component. The 
magnetic field inclination is available for both components in 
the LOS reference frame or relative to the local surface normal; 
I will use the inclination to the surface normal in the following. 
For the background field, I then determined the magnetic flux as 

®b 8 (r) = (B bg • cos y bg • f bg )(r) • A res , (3) 

where A res = (0.36 • 725) 2 km 2 is the area corresponding to a 
pixel, y bg the magnetic field inclination relative to the surface 
normal, and r measures the distance to the spot center. The fill 
factor f bg was set to 1 , since the background field is present in 
both inversion components with identical properties. 

For the flow channel component, I first converted the fill fac- 
tor to an effective size d e /f in the x-y-plane. Since the inversion 



8 



C. Beck: An uncombed inversion of multi- wavelength observations 



cannot provide information about the spatial organization of the 
two components but only the total fill factor inside the pixel, the 
definition of d e ff is ambiguous. I used two different solutions for 
d e ff. In the first, I defined d e ff\ as the edge length of a square 
whose area is equal to the fill factor by 



LOS 



2 effl 



— -yj^res ' ffc 



(4) 



like sketched in the upper right graph of Fig. [T21 This solution 
makes use of the fact that the two spatial dimensions (x, y) of 
the pixel have identical sizes of (7/36, thus a square shape is sug- 
gested. Taking into account the inclination of the fc component 
close to 90 deg (either relative to the LOS or to the surface nor- 
mal), this solution, however, provides a inconsistent topology 
since the nearly horizontal field lines of the fc component would 
have to terminate abruptly at some spatial location. A second 
solution is to demand that one of the axes (x, y) of the area that 
corresponds to the fill factor has to extend along the full size of 
the 3-D volume {bottom graph). Then the effective size in the 
plane perpendicular to the orientation of the magnetic field lines 
is given by 



2 eff2 



•ffc 



261km 



(5) 



Along the LOS, the inversion provided the stratification of 
the magnetic field with optical depth. The conversion to geo- 
metrical height in the previous section then yielded the FWHM 
of the Gaussian perturbation in km. I integrated the resulting 
Gaussian function 



Bf C (r,z) = exp - 



(z-zoY 
~2-cr 2 (r) 



■B fc (r) 



(6) 



over a height range of 1000 km in z, to cover the full extent of 
the fc component. The width <x(r) was once derived from the 
FWHM values as shown in Fig.Qj] corrected for the inclination 
of the fc component to the LOS, and once set to correspond to a 
fixed FWHM of 250 km. 

The total flux Oy c of the fc components was then derived 
from 



® fc (r) = d ef f(r) 



M00( 
JOkm 



1000km 



B fc (r,z)dz 



(7) 



for three cases, using the FWHM from the inversion and the 
two different definitions of d e ff as given above, and once with 
a fixed FWHM and d efn . The correction of the FWHM with 
the LOS inclination should provide that the effective area used 
is perpendicular to the field lines of the fc component. For the 
calculations of magnetic flux, I ignored the off-center position 
of the sunspot at about 30 deg heliocentric angle. The assumed 
size A res as given by slit width and spatial sampling along the slit 
actually corresponds to a slightly larger area (l/cos30° = 1.15) 
on the solar surface, but the projection effects enters both in the 
bg and fc flux and should thus be negligible. 

Figure [13] then displays the magnetic fluxes derived from 
the different approaches (upper panel), and the ratio between 
the flux of the bg and fc component (lower panel). Using the 
definition of d e ffi yields a flux Oj c of about 3.8xl0 17 Mx per 
pixel that is almost constant throughout the spot, whereas d e ff\ 
leads to a radial decrease of ®/ c . The flux of the bg compo- 
nent decreases with radius, due to the radial decrease of field 
strength and the increase of field inclination (see Fig. [T4l below). 
Assuming that the the use of d e ffi provides the more consis- 
tent value for Oj c , the ratio 0^/0/ c changes from about 4 at the 
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Fig. 12. Determination of the flow channels' area from the inver- 
sion results. Left top: definition of fill factor /. Top right, flow 
channel at 90 deg to the LOS. The dashed vertical line denotes 
the FWHM determined by the code. d e ff\ denotes the effective 
size of the fc component. The parallel inclined arrows denote 
the orientation of the bg field. Middle row. same for LOS incli- 
nations of the fc component of 60 and 120 deg. Bottom: another 
possible solution for the effective size d e ffi. 
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Fig. 13. Magnetic flux in the two inversion components. Top: to- 
tal flux contained in the background field and flow channel (fc) 
component. Bottom: ratio of background field to flow channel 
component flux. The horizontal line marks a ratio of unity. 



umbral-penumbral boundary to 1 at the outer white-light bound- 
ary of the sunspot. This would imply that 20% to 50% of the 
magnetic flux in a sunspot participate in the Evershed flow and 
thus presumably also in the vertical energy transport throughout 
the penumbra. 
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5. Discussion 

I applied the uncombed "Gaussian" inversion with its capability 
to reproduce the NCP to the full penumbra of a sunspot, with the 
information of five spectral lines as input. The initial model of 
the fit was taken partly from a 2C inversion with constant mag- 
netic field parameters to improve the convergence (AB, Av, AO, 
etc.), but no a priori information on the vertical structure was 
provided. The resulting best-fit spectra reproduce the observed 
NCP satisfactorily on the limb side for the VIS lines near 630 nm 
and the Fei line at 1564.8 nm, but fail partly on the center side. 
The NCP of the near-IR lines is also generally reproduced worse 
than that of the VIS lines, mainly in the amplitude of the best-fit 
NCP values that fall short of the observed values. There are some 
possible reasons for that, with the most important of course that 
in the inversion process the total least-squares deviation between 
observed and synthetic profiles is minimized, not the difference 
of the NCP, i.e., the mismatch in the NCP can be insignificant for 
the best-fit solution. The NCP of the near-IR lines also depends 
stronger on the exact location and the width of the Gaussian per- 
turbat ion since their formati on height is smaller than for the VIS 
lines (iCabrera Solana et all 12005). The Gaussian shape of the 
perturbation sets an upper limit on the steepness of the gradi- 
ents in velocity and magnetic field; i t does not correspo nd to a 
sharp discontinuity as for instance in lBorrero et al.l (120071) . The 
visible lines with their larger formation height are less sensitive 
to the exact location of the perturbation, and the gradients along 
the LOS contribute over a larger range in optical depth. The in- 
version code often used a very broad Gaussian perturbation on 
the inner center- side penumbra which then reduced the resulting 
NCP of the near-IR lines. The sign change of the NCP of the 
VIS lines in the radial direction on the center side is not repro- 
duced by the inversion, but a ring of positive (negative) NCP just 
around the umbra in the 630.37 nm line (1565.2 nm line) is. 

T07 suggested that the sign change of the NCP can be re- 
produced if the field strength in the flow channels is stronger 
by about 0.5 kG than in the background field. This proposed 
difference is in some conflict with the inversion results where 
the bg field is found to be stronger all throughout the penum- 
bra, and of equal streng t h at the outer penumbra! boundary 
dBellot Rubio et alJ 120041 iBorrero et al.l 120041: iBeckl 120081 or 
Fig. [H]). I suggest a different reason for the sign change, which 
also could explain why the Gaussian inversion was unable to 
reproduce even only the correct sign of the NCP in the cen- 
ter side. The bg component used parameters constant with op- 
tical depth, and thus, it did not contribute any NCP to the 
best-fit spectra. Without showing flow velocities as high as for 
the "flow channels", the bg component also has a significant 
velocity component of up to 2 kms" 1 in the outer penumbra 
(lowe r panel of Fig. [141 iBellot Rubio et al.|[2QQ4t IBorrero et all 
2006). The horizontal and vertical velocities have here been 
computed using the azimuthal variation of the LOS velocity all 
aro und the spot, also including the center side penumb ra (see 
e.g. ISchlichenmaier & Schmidtll200Qfc lBeckll2006L 120081) . Line- 
of- sight gradients in the flow velocity and additionally in the 
field strength of the bg component could create a non-zero NCP 
contribution from this component as well. Since the bg field lines 
are closer to being parallel to the LOS on the center side than on 
the limb side, this contribution would be more prominent in the 
Stokes V profiles of the center side and could possibly produce 
the sign change of the NCP. 

The conversion of the location of the Gaussian perturbation 
to a geometrical height scale was done with simplifying assump- 
tions on the atmospheric density stratification. A common height 
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Fig. 14. Radial variation of magnetic field strength and LOS ve- 
locity. Top: radial variation of the field strength in flow channel 
(thick grey) and background field (thin black). The difference be- 
tween them is plotted in the lower half. Bottom: horizontal and 
vertical velocities for background field and flow channel com- 
ponent from the fit of a sinusoidal to the azimuthal variation 
of vlos- Black: background field, red: flow channel component. 
Dashed: vertical velocity, solid: horizontal velocity. 



scale as used in San chez Almeidal (120051) would be an option. I 
think, however, that a revised or improved conversion between 
optical depth and geometrical height will not change the picture 
of the flow channels' topology significantly. Changes of the lo- 
cation by for instance up to ±100 km would, e.g. not remove the 
flow channel from being close to the z=0 km or change its radial 
behavior strongly. The flow channel topology from the Gau ssian 
inver sion is in good agreement with the one derived in Beck 
(2008) from the integration of the field inclination, even if the 
two methods are fully independent of each other: the location 
of the perturbation is determined separately for each pixel in the 
inversion, whereas the integration uses the radial variation of the 
inclination to derive the geometry. 

The width of the flow channel as given by the inversion code 
is more critical. On several pixels, especially in the inner penum- 
bra, the LOS cannot have penetrated through the Gaussian per- 
turbation since the lower boundary is significantly below log 
r=0. The information on the vertical extent is missing in this 
cases, the code has only used the location of the upper boundary 
to produce the gradients needed to fit the spectral lines. In these 
location, the inversion results only provide information on the 
upper limit, a maximum height of about 200 km above the z=0 
km level. 

With the same caveat that the "width" of the fc structure 
seems an ill-defined quantity in some cases, the derived mag- 
netic flux of the two components indicates that in the mid to 
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outer penumbra about 20-50% of the magnetic flux appears in 
the form of flow channels. This large fraction would imply that 
the temporal evolution of the penumbral fields has to be taken as 
the permanent and global re- arrangement of all magnetic field 
lines rather than isolated events of ascending flow channels. 

Throughout the analysis of the data presented here I now 
have used the terminology of a "background field" and a "flow 
channel" component, suggestive of an interpretation in terms of 
horizontal flux tubes embedded in a less inclined field. If this 
picture is representative of the penumbral magnetic fields is, 
however, an open question. At first, this terminology has to be 
taken simply as a convenient way of describing the two dis- 
tinct magnetic components required to reproduce the observa- 
tions, where one of them shows a larger flow velocity, a weaker 
field strength, and a larger inclination to the local surface normal 
than the other, similar to the "minor" and "major" component in 
ISanchez Almeidal (120051) . The two component s can be combined 
into a single one wit h strong LOS gradients (Math ew" et al.ll2 003: 
iBorrero et al.l l2Q04h that still reproduces the observed spectra. 
This approach runs into the problem that the magnetic field lines 
of sunspots cannot extend to the upper solar atmosphere and the 
corona since one of the components requires nearly horizontal 
fields, which is at odds with coronal observations. The anal- 
ysis results thus require two distinct magnetic components in 
the penumbra, but they do not provide directly the possibility to 
choose between any of the three penumbral models (MISMA; 
field-free gaps; flux tubes). One argument in favor of the flux 
tube picture is the spatial radial coherence of penumbral fila- 
ments over several thousand km. In the MISMA picture, it is 
hard to envisage how such a large-scale structure can be formed 
by magnetic fields s tructured on the smal lest scales without a 
common orientation. llchimoto et al.l (120071) recently reported the 
presence of roundish patches in a sunspot on disc center that ap- 
peared preferentially in the mid and outer penumbra and cor- 
respo nded to strong downflows. ISanchez Almeida & Ichimotol 
(2009) interpreted this patches in the MISMA picture, but they 
also fit well to the flux tube model. In the outer penumbra, the 
inclination of the flow chan nel compon ent shows on average a 
downward orientation (e.g., [Beck 2008), which implies down- 
flows for field- aligned mass motions. The rev isi on of the original 
simulations of Schlic henmaier et all (1 19981) in ISchlichenmaierl 
( 2002) with the peculiar "sea-serpent" shape provides additional 
indications that also in the flux tube model downflows can be ex- 
pected in the inner penumbra. The vertical velocity of the flow 
channel component in Fig. 1741 shows the same: v vert < ms" 1 
(= v is oriented downwards in this case) for r/r spot > 0.9. A 
downward oriented flux tube as drawn in the middle left panel 
of Fig. [T2l naturally would produce a roundish downflow patch 
in a horizontal cut. One caveat for the field-free gap model is 
also actually the failure of the Gaussian inversion to reproduce 
the sign change of the NCP on the limb side and the generally 
worse fit to the near-IR lines: the presence of gradients in LOS 
velocity and field strength is not suffic ient for reproducing the 
observations (IScharmer & Spruitll2006l) . only if they happen to 
be the correct gradients. 

6. Conclusions 

The uncombed "Gaussian" inversion is able to reproduce the 
NCP of simultaneous observations in VIS and near-IR spec- 
tral lines well on the limb side of a sunspot, where the signa- 
ture of the two different magnetic components in the penum- 
bra is clearest. Even with a mismatch in the NCP, it still is 
able to reproduce the observed spectra of five spectral lines 



satisfactorily all throughout the penumbra, with at least the 
same quality as a 2-component inversion with constant mag- 
netic field properties (Figs. IB. II to IB. 3b . The inversion setup 
can be interprete d as embedded flow channe ls like in the pic- 
ture suggested bvlSolanki & Montavonl (1 19931) . As remarked by 
IScharmer & S pruit (20061). the agreement does not prove the cor- 
rectness of the flux tube model, but it still proves that the model 
definitely is not at odds with the observations. 

Around 20-50% of the total magnetic flux in the penum- 
bra shows the characteristics of the component addressed as 
flow channel, implying a permanent re-organisation of all mag- 
netic field lines in the penumbra. This could maybe allow to re- 
plenish the ra diative losses of the penumbr a also in the moving 
tube model of ISch lichenmai er et al.l (Il998l) . The calculations of 
ISchlichenmaier et al. ( 1999]) yielded an energy supply by a sin- 
gle flux tube that was insufficient to compensate the penumbral 
energy losses on a radial cut, but eventually the rate and/or num- 
ber of such tubes should be increased until their magnetic flux 
reaches 50% of the static background field component. 

Two peculiarities of the NCP observed in the sunspot at 
30° heliocentric angle stand out prominently: a sign change of 
the NCP of the VIS lines on the center side as reported by 
iTritschler et a l. ( 2007), and a ring of positive and negative NCP 
just around the umbra in the Tii line at 630.37 nm and in the 
Fei line at 1565.2 nm, respectively. The first item is not repro- 
duced by the inversion and is also missing in the theo retical 
calculations of, for insta nce, iMtiller et aD (120021 12006b . or in 
ISanchez Almeidal (|2005[ Fig. 19, for a sunspot closer to disc 
center). It would be interesting to see i f this sign change is re- 
produced by the flux tube model used in lBorrero et al.l (120071) or 
Bor rero & Solankil (1201 Ol) on radial instead of azimuthal paths. 
The ring pattern around the umbra might be well suited for 
comparisons between observations and the MHD simulations of 
iRempel et al.l (l2009allbh . since the penumbral "filaments" in the 
simulations are somewhat shorter than any observed filaments, 
but provide information on the inner footpoints near the umbral- 
penumbral boundary. 

The observations used in the present study were taken at the 
VTT in 2003, when it was only possible to improve the spatial 
resolution by a correlation tracker providing image stabilization, 
and with a 50% loss of light due to the use of an achromatic 
beamsplitter to feed the two VIS and near-IR instruments. With 
the growing activity of the new solar cycle, it would now be 
possible again to obtain improved multi- wavelength data sets at 
the VTT for future studies, with both a higher signal-to-noise 
ratio thanks to a new dichroic BS and better spatial resolution 
thanks to adaptive optics. 
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Appendix A: Changes between initial and best-fit 
model atmosphere 

The Gaussian inversion was initialized with the difference be- 
tween the two inversion components in the 2C inversion as am- 
plitude of the Gaussian perturbation. Figure IA.1I shows how 
much the initial value was modified in the inversion process. The 
changes in the difference of the LOS magnetic field azimuth, 
AO, the difference of the LOS magnetic field inclination, Ay, 
or the difference in field strength, AB, were actually minor. The 
field orientation (AO, Ay) changed more than the field strength, 
but this could also be due to the fact that the values for the sec- 
ond component with the Gaussian perturbation were taken at a 
fixed optical depth of log r = 0. 

Appendix B: Profile examples 

Figures IB. II to IB. 31 show the spectra of six different locations 
inside the penumbra, marked by crosses in Fig. [T] The best-fit 
profiles of the 2C inversion (blue lines) can be seen to deviate 
strongest from the observations in the case of Stokes V of the 
VIS lines like for instance in the left panel of Fig. IB. 21 (profile 
no. 3), while still well reproducing the near-IR V spectra at the 
same time. 

The best-fit atmosphere models in Fig. lB.4l show that in some 
cases (profiles no. 4 and 5) the Gaussian perturbation actually is 
converted to a shape quite different to a Gaussian, i.e., an ex- 
tremely broad Gaussian where the lower boundary of the per- 
turbation is located far below logr = 0. In these cases, the 2C 
inversion setup with constant field parameters and the Gaussian 
inversion are actually as good as identical approaches to analyze 
the spectra, and hence yield very similar best-fit spectra. 
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Fig. B.l. Profile examples showing the observed spectra {black crosses), the best-fit of the 2C inversion (blue line), and the best-fit 
of the Gaussian inversion (red line). In each panel, Stokes IQUV are shown from left to right, and the lines 1564.8 nm, 1565.2 nm, 
630.15 nm, 630.25 nm, and 630.37 nm from bottom to top. The locations of the profiles inside the FOV are marked by crosses in 
Fig. HI 
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Fig.B.4. Atmospheric stratifications of the best-fit result of the Gaussian inversion for the profiles shown in Figs. IB.ll to IB.3I The 
layout of each panel is identical to that of Fig. [2] 
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